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Abstract 

> 

[--~. ' Due to the failure of Birkhoff 's theorem, black holes in /(i?) gravity theories in 

CTN ■ which an effective time-varying cosmological "constant" is present are, in general, 

dynamical. Clifton's exact spherical solution of B}^ gravity, which is dynamical 
^ and describes a central object embedded in a spatially flat universe, is studied. 

f~^ ' It is shown that apparent black hole horizons disappear and a naked singularity 

^D . emerges at late times. 
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1 Introduction 

The study of Type la supernovae [l3lliTllllll2lll5l[6Tl[30lll6l[3] revealed that the 
universe is currently accelerating its expansion. This discovery has generated an enor- 
mous amount of activity and theoretical models in order to find an explanation of this 
phenomenon. The most common models are based on General Relativity (GR) and 
invoke mysterious forms of dark energy (see [32] for a list of references). However, dark 
energy, possibly even phantom energy, is too exotic and ad hoc and attempts have been 
made to model the cosmic acceleration without dark energy. f{R) theories of gravity 
reminiscent of the quadratic corrections to the Einstein-Hilbert action introduced by 
renormalization have been re- introduced in the metric dU] , Palatini [53] , and metric- 
affine [HH [321 IS31 [SSI [SZ] formulations and have received much attention in recent years 
(see [SSI [16] for reviews and [SHI [SH ISl [201 [3 ISHl \M\ for introductions). Emilio Elizalde 
has given many contributions to the development of f{R) gravity and cosmology. 

While many cosmological and other aspects of f{R) gravity (stability, weak- field 
limit, ghost content) have been discussed in recent years, it is important to understand 
spherical solutions (both vacuum and interior) in these theories [2S1 HI [SOI [IS [SB [21 HSl 
[^[62]. 

Metric f{R) gravity is described by the action 

S = ^ f d^'x^/^fiR) + S^'""**"") , (1.1) 

2k J 

where f{R) is a non-linear function of its argument and S^"^""*^^^^ is the matter part of 
the action. R denotes the Ricci scalar of the metric Qab with determinant g, k = SnG 
where G is Newton's constant, and we adopt the notations of Ref. |64] . 

The Jebsen-Birkhoff theorem of GR fails in these theories, adding to the variety of 
spherical solutions [23]. Of particular interest are black holes in generalized gravity, 
which have been studied especially in relation to their thermodynamics ([21] and ref- 
erences therein). Since f{R) theories are designed to produce an effective dynamical 
cosmological constant, physically relevant spherically symmetric and black hole solu- 
tions are likely to describe central objects embedded in cosmological backgrounds. This 
kind of solution is poorly understood even in GR, although a few examples are avail- 
able there [SSl[2Il[Sl[Hll3Sl[3Sl[2Sl[lZl[2Sllini[3a[2S]. Even less is known about f{R) 
black holes, which are certainly worth exploring. Here we consider a specific solution 
of vacuum f{R) = R^^^ gravity discovered in [11]. Solar System experiments set the 
limit^ (5 = (-1.1 ± 1.2) • 10"^ on the parameter 5 [IH [H [121 [TSl [6S] , while local stability 
requires f"{R) > [TTl [TH [TH [38] , hence we restrict to the range Q < 5 < 10"^. 



^See also [14] for this specific form of the function f{R). 



The solution of [TT] is dynamical and describes a time-varying central object embed- 
ded in a spatially flat universe in vacuum R^^^ gravity. This solution is made possible 
by the fact that the fourth order field equations of metric f{R) gravity 

f'{R)Rat -^9ab = V,V J'(i?) - gabOfiR) (1.2) 

in vacuo can be rewritten as effective Einstein equations 

(f- RfMl 

(1.3) 



Rab — 77 9abR 



VaV.f - gabOf + gab^^-^^ 



2 ^"" f'{R) 

with geometric terms acting as effective matter on the right hand side. This time- 
varying effective matter invalidates the Jebsen-Birkhoff theorem and can propel the 
acceleration of the universe. An equivalent representation of metric f{R) gravity as 
an a; = Brans-Dicke theory with a special scalar field potential reveals explicitly the 
presence of a massive scalar degree of freedom f'{R) responsible for these effects [55] . 
Since analytical spherical and dynamical solutions of f{R) gravity in asymptotically 
Friedmann-Lemaitre-Robertson- Walker (FLRW) backgrounds are harder to find than in 
GR (where only a few are known anyway), Clifton's solution is particularly valuable. 

2 Clifton's solution and its horizons 

In this section we describe the Clifton solution [TT] and the work [22] locating the horizons 
of this solution. 

The spherically symmetric and time-dependent solution of vacuum R^'^^ gravity of 
[TT] is given by 



ds^ = -A2ir)dt^ + a\t)B2ir)[dr^ + r^d9^ + sin'9dip^)] , (2.1) 

and (using the isotropic radius and the notations of [11\) 



B^ir) = (l + — ) A2{r)'^+''~\ (2.3) 



C2 



r 



S(1+2S) 

a{t) = t^=^ , (2.4) 



q' = l-26 + A6\ (2.5) 



Once S is fixed, two classes of solutions exist, corresponding to the sign of C2qr. The 
line element (12.1 p reduces to the FLRW one if C2 — )■ 0. In the limit 5 — )■ in which 
the theory reduces to GR, eq. (12. ip reduces to the Schwarzschild metric in isotropic 
coordinates provided that C2qr > 0, hence positive and negative values of r are possible 
according to the sign of C2, but we assume r > 0, C2 > and take the positive root in 
the expression q = ±y/l -26 + 46"^, so that g ~ 1 - 5 as 5 ^ 0. The solution f l2T^ -f E3]l 
is conformal to the Fonarev solution |27] which is conformally static [33] , and therefore is 
also conformally static, similar to the Sultana-Dyer [591 EB M E] and certain generalized 
McVittie solutions [26] of GR. 

In order to identify possible apparent horizons, it is convenient to cast the metric (12. ip 
in the Nolan gauge. Using first the Schwarzschild- like radius 
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r = r I 1 + 
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(2.6) 



giving dr 



1 - % 



-1 



dr and then the areal radius 
a{t)^B2{r)f 



P = 



a{t) f A2{r 



2 



(2,7) 



(1 + ^)^ 
the line element (12. ip takes the form 

ds"" = -A^dt^ + a^Af~^df^ + p'd^^ . (2.8) 

Denoting the differentiation with respect to time with an overdot and using the identities 



df 



dp — A2 



9+2(5-1 



df dt 



dp — A2 



q+2S~l 



dr dt 
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0+2(5-1 , , , 25-1-q 
/(" 2 1 2{q+25~l)C2 4 2 
2 ~^ q f 2 



9+2(5-1 



aA^ - 



-C{r) 



(2.9) 



one obtains 
C(r) 



^ ^ 2{q + 25-l) C2 ^^.q 
q r ^ 



1 + 



2(g + 26-1) C2a 



q P 

which turns the metric into the Painleve-GuUstrand-like form 

,2(5-1) 



2(5-l-(j 
A2 ' 



(2.10) 
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-A, 



dp^ 
AlC^ 



1- 



A 



C2 



• 2-2 

a r 



dt'- 



-q+26-l 

2Ao ' 



^2 



a r dtdp 



+ p'dn\ 



(2.11) 



Now we introduce a new time coordinate t defined by 

1 



dt 



F{t,p) 



[dt + /3{t,p)dp] 



(2.12) 



in order to eliminate the cross-term dtdp. Here F{t,p) is an integrating factor wliicli 
guarantees tliat di is an exact differential and is determined by 



dp \FJ dt \F 



The line element becomes 
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The choice 



/3 



-q+2<5-3 
^2 ' 
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n2 42(^-1) 

<-^ 1 _ ^ A2^2 



C2 
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eliminates the dtdp term and casts the metric in the Nolan gauge 



ds^ = -A2DF^dP + 



AlC^ 



<?-l u2 J2 



1 + 



A,'-'H'' p 



dp^ 



+ p^dn'^ , 

where H = a/ a is the Hubble parameter of the background universe and 



D = 1 
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2(5-1) 
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a r = 1 
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(2.13) 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



Using the second of these equations, the hne element (12.1 6p assumes the simple form 



ds' 



-A2DF'dt' + 



AlC^D 



+ p'd^' 



(2.18) 



The apparent horizons, if they exist, are located at g'''' = 0, which yields A^C^D = 
and Al (C2 _ i72^2^-9-i) = 0. Therefore, g^^ vanishes if A2 = or H'^R'^ = C^Al+\ 
A2 vanishes at r = C2, which describes the Schwarzschild horizon when 5 — )■ (the 
GR limit). This locus corresponds to a spacetime singularity because the Ricci scalar 

R = y^ (^.) diverges as r — )■ C2 (it reduces to the usual FLRW value 6 (H + 2H'^ j 
as C2 — !■ 0). This singularity is strong according to Tipler's classification [BO] because 

q + 2S-l 

the areal radius p = arA2 ^ vanishes when r = C2 for 5 > 0, in contrast with the 
Schwarzschild metric corresponding to 5 = in which p = r = 46*2 at r = C2. 
The second possibility H^p"^ = C'^A'^ yields 



Hp=± 



1 + 



2(g + 26-1) C2a 



a: 



25-1-q 
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2+1 



(2.19) 



with the positive sign corresponding to an expanding universe. When 5 — )■ 0, this 



equation reduces to Hp 



1 + 



2SC2a 



A, 



2 y 



Al-'. 



To gain some insight, consider the following two limits. As C2 — t- (the central object 
disappears and the solution is FLRW space), r = r and p become a comoving and a 
proper radius, respectively, while eq. (I2.19P reduces to Hp = 1 with solution pc = l/H, 
the radius of the cosmological horizon. In the limit 5 — )■ in which the theory reduces 
to GR, eq. (|2:T9|) reduces to A2 = or r = ^3 with H = 0. 

Using eqs. (12.41) and (12. 7p . the left hand side of eq. (I2.19P is expressed as 
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(2.20) 



where x = C2/T, while the right hand side of (I2.19P is 

2{q + 2d-l) X 



3±1 r- 

1 — a; \ 1 



1 + x 



1 + 
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(2.21) 



Eq. fl2.19p then becomes 
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(2.22) 



(note that ^-fj-/*^" is positive for < 5 < ^|^ ~ 0.366). 

At late times t, the left hand side of eq. (I2.22p vanishes, x ~ 0, and there exists 
a unique root of the equation locating the apparent horizons, which corresponds to a 
cosmological horizon, consistently with the fact that r — )■ cxd as x = C2/r — )■ 0. The 
limit a; — )■ can also be obtained when the parameter C2 — ?■ , in which case Hp — )■ 1 
and r ^ p ^ H~^ = g,|~ ^. t is the radius of the cosmological horizon of the FLRW 
space without a central object. Hence, there is only a cosmological apparent horizon and 
no black hole apparent horizons at late times: the central singularity at p = becomes 
naked. 

The radii p of the apparent horizons and the time t can be expressed in the parametric 
form 



p{x) 



t{x) 



t{x) 
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(2.23) 



,(2.24) 



using X as a parameter. Fig. 1 reports p versus t for the parameter values C2 = 1 
and 6 = 0.13, showing that two inner horizons develop after the Big Bang covering the 
central singularity p = 0, then they approach each other, merge, and disappear, while a 
third, cosmological horizon keeps expanding. The p = singularity becomes naked after 
this merging event. 



3 Discussion and conclusions 

Cosmologists may be detecting deviations from GR and therefore it is necessary to 
understand spherical solutions of f{R) gravity, which has been proposed as a simple 
alternative to the mysterious dark energy. Since the Jebsen-Birkhoff theorem fails in 
these theories, spherical solutions do not have to be static. f{R) theories are designed 



with a built-in dynamical cosmological constant to model the present acceleration of 
the universe, hence analytical spherical solutions describing a central object embedded 
in a FLRW background are the relevant ones. Unfortunately, such solutions are poorly 
understood even inGR[59l|2ll|9l[8l|35l|36l[25lll71[29lll0llia|26]. It seems difficult to 
find generic solutions describing black holes embedded in FLRW backgrounds. Finding 
numerically spherical interior solutions of f{R) gravity is also an active area of research 
[SSlKEnillSlEIlElIISlEZlES]. All these issues deserve further attention in the future. 
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Figure 1: Radii of the apparent horizons of Chfton's solution (vertical axis) versus time 
(horizontal axis) for the parameter values C2 = 1 and 5 = 0.13. 



12 



